Stability Condition of a Strongly Interacting Boson-Fermion Mixture across an 

Inter-Species Feshbach Resonance 
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We study the properties of dilute bosons immersed in a single component Fermi sea across a broad 
boson-fermion Feshbach resonance. The stability of the mixture requires that the bare interaction 
between bosons exceeds a critical value, which is a universal function of the boson-fermion scatter- 
ing length, and exhibits a maximum in the unitary region. We calculate the quantum depletion, 
momentum distribution and the boson contact parameter across the resonance. The transition from 
condensate to molecular Fermi gas is also discussed. 



Feshbach resonances are an important tool for achiev- 
ing strong interactions in degenerate quantum gases. The 
many-body system at resonance, if stable, is of great in- 
terest since the s-wave scattering length a s diverges and 
the system cannot be treated by conventional perturba- 
tive means. On the other hand, universal properties arise 
at resonance and greatly simplify the problem. However, 
in order for a degenerate atomic gas to be stable, two 
conditions are required. First, the atom loss rate should 
be small to ensure long enough lifetime and second, the 
system should have a positive compressibility against me- 
chanical collapse. 

So far, two-component Fermi gases are the only sys- 
tems whose stability at resonance has been firmly estab- 
lished and their properties have been studied extensively 
in the last few years In contrast, bosons at resonance 
suffer from rapid atom loss [2| and collapse instability [3j , 
and are not stable. 

For a boson-fermion mixture, the atom loss rate (domi- 
nated by three-body recombination) depends on the con- 
centrations of the two species. In typical situations 
the boson density Ub is much larger than the fermion 
density n F and thus atom loss is also significant at res- 
onance. However, in the limit of low boson concen- 
tration, most three-body processes involve two identi- 
cal fcrmions and one boson, and the loss rate can be 
greatly suppressed. As for stability against collapse, 
the weak coupling mean-field theory for a uniform mix- 
ture gives kpa-BB ? (&FaBF) 2 (l + l) 2 /{^l) @|) where 
k F = (67t 2 ?if) 1 / 3 , 7 = me/mp, obb and obf are the scat- 
tering lengths between bosons, and between bosons and 
fermions, respectively. In addition to the 6 Li- 23 Na and 
40 K- 87 Rb mixtures studied previously @, new mixtures 
have been very recently been realized in the laboratory, 
including 6 Li- 87 / 85 Rb7] and 40 K- 41 K [8]. It is therefore 
very important to generalize the stability condition of 
boson-fermion mixtures to the strongly interacting region 
where obf ±oo. 

There have been a number of theoretical studies of 
boson-fermion mixtures. Most address the region away 
from unitarity 0, [13] or deal with resonance physics us- 



ing a mean- field treatment of a two-channel model 11]. 
Some recent works studied boson-fermion pairing effects 
within a single-channel model for a broad resonance 12 1 . 
However, the interaction between bosons, which is cru- 
cial for the stability of the mixture as we will show, is 
ignored in most of these analyses, and thus the question 
of stability is still an open one. 

The Hamiltonian % for a boson-fermion mixture can 
be separated into Hq and ri\ (h = 1), 
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where Ubf and Ubb are the zero-range pseudo-potentials 
that produce the scattering lengths obf and clbb ■ In 
this work we consider the situation where cjbf can be 
tuned by a broad Feshbach resonance and the boson- 
fermion interaction energy can be of the same order as 
the Fermi energy, and hence cannot be treated pertur- 
batively. In this work we will use a Jastrow-Slater vari- 
ational wave-function and the lowest order constrained 
variational (LOCV) approximation to study the ground 
state properties of Hq. On the other hand, we shall re- 
strict ourselves in the regime x = iib/hf <C 1 so that 

1/3 

n B a,BB is always small and positive. Therefore we can 
treat ri\ as a perturbation and evaluate it at the mean- 
field level. Two main results of this Rapid Communica- 
tion are summarized as follows: 

(i) The energy density of the mixture is found to be 
given by 
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-^A( V )nBe F + ^gBBn 2 B [l + D{ V )], (2) 



where ep = fc F /(2mp), #bb = 47tcibb / tub and r\ = 
I/^fObf). A(r/) and D(rf) are two universal functions of 
r] and are shown in Fig. Q] These two functions have sim- 
ple physical interpretations. In the boson-fermion mix- 
ture, minority bosons are dressed by a polarized cloud of 
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FIG. 1: (Color online) Solid line: Universal function A{rf) and 
D(n) as a function of coupling constant n = I/^fosbf). The 
dashed and dash-dotted lines in (a) are the mean-field value 
and molecular binding energy in the two limits, respectively. 
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FIG. 2: (Color online) (a) Momentum distribution of non- 
condensed bosons A r B(k)fc 2 /fcp for different n — I/^fosbf); 
(b) Contact C (in unit of fcp 4 ) as a function of 77. For both 
cases we set x — bb/jif = 0.1. 



fermions and form bosonic polarons; A{rj)t-p is the po- 
laron binding energy for 7=1. D(j]) is the fraction of 
non-condensed bosons (quantum depletion) arising from 
boson-fermion interactions. 

(ii) With the ground state energy in Eq. ([2]) we deter- 
mine the mechanical stability condition, 



^FOBB > Cc(??,7)^)) 



(3) 



across a Feshbach resonance. In the dilute limit x <C 1 
to zeroth order in x, 
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Here £ c is a universal function of 77 for given 7 as shown 
in Fig. (Ha). It attains a value of the order unity at 
resonance and is maximum in the unitarity region. 

J astrow- Slater Calculation of Energy Density. We use 
a variational wave-function 
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where |$fs) represents a Fermi sea with iVp fermions 
and /(r) is the Jastrow function describing the two- 
body correlations between bosons and fermions. /(r) 
deviates from unity only within the so-called healing 
length d, i.e., f(r > d) = 1. To estimate the en- 
ergy £ = (^jslHol^jsl/i^jsl^Js), we use the so-called 
LOCV approximation [13[, which basically treats corre- 
lations to lowest order (0(f 2 — 1)) in a linked cluster 
expansion. LOCV was first used to study 4 He |13j . and 
has also recently been applied to study bosons with large 
scattering lengths [l4| and Fermi gases at resonance [l5| . 
The calculation is both simpler and more transparent 
than Monte- Carlo simulations while providing a fair ap- 
proximation to Monte-Carlo results. 

Within LOCV, we find E /V = +£ B f, where £{? = 
3np£F/5 is the energy of the free fermions and £ B f is the 
interaction energy 



£bf = n B n F / d J r /(r) 
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between bosons and fermions with m r = vti b ttiy / {jn B + 
rnp). Eq. ([6]) can be rewritten as £bf = An B = 
A(ri, "f)n B e F , where A is found by solving [ - 5^73^ + 
U BF (r)] rf{r) = Xrf(r), subject to the constraint 
A-ku f Jq dr r 2 \f(r)\ 2 = 1 an d the boundary conditions 
f(d) = 1 and f(r)'\ r=d = 0. It fol- 

lows that the universal function A is given by A(n, 7) = 
HV:l)/ e F = (1 + l)A(n)/(2j). For a given value of 77, 
there is only a change of prefactor for different mass ra- 
tios 7. 

The numerical solution of A(rj) across the resonance is 
shown in Fig. [lja). In the limit r\ — > —00, we find A(r\) = 
4/(3nri), so the interaction energy density reduces to the 
mean-field result 2tt a BF n F n B /m r . In the opposite limit 
i] — > +00, A{rf) = —2n 2 , and A becomes the binding 
energy of the molecule — l/(2m r a BF ). 

At resonance (77 — 0), we obtain A = —0.64 for 
equal masses. It is interesting to compare this result 
with related studies of highly polarized fermion-fermion 
mixtures, where minority fermions form polarons and 
the system behaves as a polaron Fermi liquid. Various 
theoretical approaches have been used to estimate the 
binding energy of a single impurity IB- 3~sj] . A varia- 
tional approach gives A = —0.61 [17l| and diagrammatic 
Monte Carlo give s A = —0.62 [18j. Experiments yield 
A = -0.64(7) lH. 

LOCV can also be used to determine the quantum de- 
pletion of bosons ndep as well as their momentum distri- 
bution N B (k): 
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The depletion fraction D(n) across the resonance is plot- 
ted in Fig. [ljb). It is a monotonically increasing function 
of rj. At resonance, we find D = 0.13. Similar calcula- 
tions have been carried out for a resonant Bose gas by 
Cowell et al. 14|. The momentum distribution of de- 
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FIG. 3: (Color online) Phase diagram of dilute bosons im- 
mersed in a single component Fermi gas: (a) in terms of 
77 = l/(fcpaBp) and £ = kpass', (b) in terms of iiFa| B and 
clbb/clbf- The dashed line in (a) is the mean- field result in 
the weak coupling limit, (c) Critical value Cc as a function of 
boson concentration x — ne/fiF for different rj. For all cases 
we set tub = mp. 



pleted bosons are plotted in Fig. [2fa) for different n, 
and the contact C = lim^+oo fc 4 A r B(k) [20| as a func- 
tion of i] is plotted in Fig. HJb). We note that at low 
momentum, A r e(k) will be modified by boson-boson in- 
teractions. However, this will not affect our analysis of 
the large-k limit of Nb (k) or the stability analysis below. 

With the solution for |\Pjs), the energy correction due 
to "Hi is obtained to lowest order: 



E x 1 (* JS |fti|*js) 1 ari , . n ,o , Q , 
= -gBBn B [l + 4D{T])\. (9) 



V V (*j S |*js) 



Physically, the first term in Eq. ([9]) is the Hartree en- 
ergy between the interacting bosons, and the second term 
can be interpreted as the Fock energy and pair hopping 
energy arising from interactions between condensed and 
non-condensed bosons. 

Stability Condition. Using Eqs. (|6) and ©, Eq. © 
gives the total energy of the system. In the weak coupling 
limit our theory agrees quantitatively with the perturba- 
tion results of Ref. [lOj . From Eq. © we obtain 
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where C = fepiBB- The stability of the system requires 
(i) djip/dnp > 0; (ii) djis/dnB > and (hi) the stability 
against small density fluctuation requires Det(M) > 0, 



M = 



djip/dnp djip/dnB 
d^LB/dnF 9/iB / Oub 
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Here (i) and (ii) are easy to satisfy, while (iii) gives 
C > (ciVili 30 )' an d Cc is given by Eq. Q in the di- 
lute limit x <C 1. For equal masses, the phase dia- 
gram in the dilute limit is shown in Fig. [3] Physically, 
collapse is a global instability due to negative boson- 
fcrmion mean-field energy, and only sufficiently strong 
boson-boson repulsions prevent the density from increas- 
ing without bound. In the 77 — ( plane (Fig. EJa)), 
the phase boundary is given by the universal function 
Cc(??)- At resonance, we find Cc — 0.43. In the limit 
7/ — > —00, one finds dA/drj = —A/j] and D(it) — > 0. 
Thus, Cc = (I+7) 2 / (27T7?7 2 ), which agrees with the mean- 
field result Cc - > [5j. In the opposite limit, rj — > +00, 
one also finds Cc — > because 2A(rj) — rjdA(rj)/drj — » 0. 
Therefore, a maximum of Cc is expected in the unitary 
region. For equal masses, the maximum of Cc is 0.49, 
above which the system is stable for all values of rj. 
Physically, the non-monotonic behavior of Cc arises from 
the off-diagonal term djiB/dnp in M. For n — > —00, 
\dfiB/dnF\ vanishes linearly with |(zbf| — > 0; while for 
rj — > +00, the attraction becomes so strong that only 
short-range physics matters, jiB eventually approaches 
the molecular binding energy — l/(2m r a BF ) independent 
of n F , and \d / dn-p\ vanishes again. Hence \d^/dn^\ 
must exhibit a maximum in between which gives rise to 
the maximum of Cc ■ 

Another notable feature of the phase diagram is that, 
for a given scattering length obb, the stability require- 
ment for 71f is opposite in the weakly interacting and 
resonance regions. For weak interactions, stability re- 
quires n-p ^ [27raBB7/(aBF(l + 7) 2 )] V^"" 2 ); while at 
resonance it requires n-p > (Cc/aBB) 3 /(67r 2 ). In Fig. [3] 
(b), the phase diagram is plotted in terms of n F a B B an< ^ 
«bb/obf- For |aBB/aBF| > 0.18, the system is stable for 
all values of n-p and for |<xbb / ^bf | < 0.18, as np increases, 
the system first becomes unstable and then becomes sta- 
ble again. In Fig. [3jc), we see that Cc depends weakly on 
the boson concentration x for small x. 

Applicable Region of Our Theory. For rj 3> 1, strong 
pairing fluctuations lead to the formation of fermionic 
molecules, and a Jastrow-Slater type wave-function is no 
longer applicable. The mean-field energy of a uniform 
atom-molecule mixture is given by [2l| 
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where omf is the molecule-fermion scattering length 
22]. For unbalanced fermion-fermion mixtures, a sim- 
ilar mean-field energy has been found to be quite ac- 
curate compared with Monte Carlo results when n > 1 
[23j | . Here we estimate the critical value rj c of the phase 
boundary by comparing the energy in Eq. pip with the 
energy Eq. @ discussed above. To zeroth order in x, n c 
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FIG. 4: (Color online) (a) SE, defined as the energy differ- 
ence between the condensate [Eq. ([2])] and atom-molecule 
mixture [Eq. energies as a function of n = l/(feFOBp)- 

We consider the case of equal masses as well as the mixtures 
of 87 Rb- 40 K and 87 Rb- 6 Li. (b) The estimated transition point 
?7 C as a function of x = iib/bf for 7 = 1 and different 



is determined by the equation 



1 + 7 2 

-Vc 



2(2 + 7) «mf 



(12) 



7 37r(l + 7)?7 C obf 

The energy comparison is shown in Fig. E{a) . For 7=1 
we obtain r/ c = 1.19, and 7y c moves toward resonance as 
7 increases. Fig. @|b) shows that ?7 C depends weakly 
on x. We note that the energy comparison only gives a 
rough estimate of the phase boundary. The nature of the 
transition, i.e. whether the loss of condensation and the 
appearance of a molecular Fermi surface occur at one 
first-order critical point or two separated second-order 
critical points, is left for future work. 

Experimental Realization. In this work, we analyze the 
stability condition and determine the phase diagram of a 
boson-fermion mixture across a Feshbach resonance. Ex- 
perimentally, at a given magnetic field near the boson- 
fermion resonance, cibb is fixed. However, if one com- 
bines other control techniques, for instance, using opti- 
cal Feshbach or microwave-induced resonances, one can 
also tune obb independently as proposed recently in Refs. 
[jij ]. The mechanical stability condition Eq. ([3]) pre- 
dicted here will provide a useful guide for experiments 
searching for stable, strongly-interacting boson-fermion 
mixtures. The depletion fraction, momentum distribu- 
tion and contact calculated here can be measured by 
various techniques such as B ragg spectroscopy [25| and 
radio-frequency spectroscopy |26| . 
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